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Abstract 

In this paper we discuss our experiences with Navier- 
Stokes (NS) codes using central differencing (CD) and 
scalar artificial dissipation (SAD). NS-CDSAD codes 
have been developed by Jameson, Martinelli, Swanson, 
and Vatsa among others. Our results confirm that for 
typical commercial transport wing and wing/body con- 
figurations flying at transonic conditions with all tur- 
bulent boundary layers, NS-CDSAD codes, when used 
with the Johnson-King turbulence model, are capable 
of computing pressure distributions in excellent agree- 
ment with experimental data. However, results are not 
as good when laminar boundary layers are present. Ex- 
haustive 2-D grid refinement studies supported by de- 
tailed analysis suggest that the numerical errors asso- 
ciated with SAD severely contaminate the solution in 
the laminar portion of the boundary layer. It is left as 
a challenge to the CFD community to find and fix the 
problems with Navier-Stokes codes and to produce a NS 
code which converges reliably and properly captures the 
laminar portion of the boundary layer on a reasonable 
grid. 


to other available methodologies. Thus JT-NS codes 
have achieved a wide following in the CFD community. 

In this paper we will present some applications of JT- 
NS codes to 3-D and 2-D problems of aerodynamic in- 
terest, including wing/body, nacelle, airfoil and multi- 
element airfoil configurations. 

We will begin our discussion with an account of 
the relative success of JT-NS codes applied to 3-D 
wing/body configurations with turbulent flow. We will 
follow with a somewhat sadder tale for 2-D airfoils in- 
volving runs of laminar flow. Our attempts to locate 
the problems with 2-D JT-NS have included detailed 
grid refinement studies which indicate numerical prob- 
lems particularly in the laminar portion of the bound- 
ary layer. These numerical problems are discussed at 
length. The inability of JT-NS codes to properly cap- 
ture the laminar portion of the boundary layer (on a 
reasonable grid) prevents us from including a stability 
analysis needed to predict the onset of transition. We 
give examples where transition prediction is very im- 
portant, including flow around a high-lift multi-element 
airfoil configuration and for around a nacelle. 


Introduction 

Boeing’s recent acquisition of a CRAY Y-MP has en- 
abled us to perform definitive grid-refinement studies 
with NS codes. We will focus attention on Jameson- 
technology (JT) codes developed by, among others, 
Jameson [1], Martinelli [2], Swanson [3], and Vatsa [4]. 
JT-NS codes employ central differences (CD) with 
scalar artificial dissipation (SAD). From the point of 
view of accuracy in a Navier-Stokes calculation, CDS AD 
is thought by some not to be as good as other available 
methods. Indeed, van Leer [5] boldly states that there 
is no hope for the flux formula of the Jameson type. 
Nonetheless, the combination of CDSAD with Runge- 
Kutta time marching, augmented with implicit residual 
smoothing and multigrid, have given JT codes a well 
deserved reputation of being fast and reliable compared 


Wing and Wing/Body Analysis and 
Design 

In this section we will compare the capabilities of the 
JT-NS code TLNS3D developed at NASA Langley [4], 
using the Johnson-King turbulence model [6] , with our 
traditional viscous/ in viscid coupled code A488. 

In figure (1) we show a comparison, for a supercrit- 
ical wing near design conditions, between experimen- 
tal data, TLNS3D, and A488. The TLNS3D solution 
matches well with test data, whereas the A488 solution 
places the shock too far back. For the many test cases, 
TLNS3D has proven to be consistently more accurate 
than A488 [7, 8]. 

Another advantage TLNS3D enjoys over A488 is the 
ability to predict flows at off-design conditions involving 



flow separation. Figure (2) shows wing upper surface 
streamlines and pressure distributions at 65% semispan 
station. TLNS3D properly predicts the trailing edge 
separation and detailed surface pressures. 

Encouraged by these successful analysis Tuns usings 
TLN53D, wort has begun on design. An iterative de- 
sign method that allows the designer to prescribe de- 
sired pressure distributions together with geometry con- 
straints, such as thickness and trailing edge closure is 
under development. Preliminary results based on this 
method are shown in figure (3). Here the target geom- 
etry is the ONERA M6 wing, and the pressure distri- 
butions are given. Beginning with a NACA 0012 wing 
section as an input geometry, the target geometry is ac- 
curately recovered within 20 design cycles. A more de- 
tailed description of the design method is under prepa- 
ration [9]. 

Practical CFD Assessment for Wing/Body 

Generally speaking, the ability of TLNS3D to prop- 
erly predict the pressures at both cruise and slightly 
off-design conditions is good. The main problems 
are laminar flow predictions and accurate drag pre- 
dictions. Accurate drag predictions are crucial to de- 
sign/optimization. Indeed, one of the design goals is to 
maximize the lifl-to-drag ratio (under constraints). The 
designer will make a considerable effort to reduce drag 
by even as little as 1%. It is estimated that a 1% drag 
reduction, for a long-range airplane such as the Boe- 
ing 777, will save the airlines 6 billion dollars, based on 
a 2,000 airplane fleet operating over a 20 year service 
life [10]. Customer airlines require that tight perfor- 
mance guarantees be offered years before the airplane is 
actually built. In this tough commercial environment, 
the accuracy and reliability requirements must be very 
high if CFD is to be depended on to help fine tune final 
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Figure 2: 3-D Navier-Stokes Streamlines and Pressures 
at 65% semispan 
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Figure 3: 3-D Navier-Stokes Design 


designs and to establish meaningful performance guar- 
antees. 
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2— D Airfoil Studies _ nR 


The most direct 2-D equivalent to TLNS3D is the 
JT-NS code FLOMGE developed by Swanson [3] which -0.2 
also incorporates the Johnson-King turbulence model. 

For some flow situations, FLOMGE gives reasonable re- 
sults. An example involving EAE 2822 case 6 is shown ct 0.2 
in figure (4). 
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Figure 5: NACA 0012, experiment (Symbols) versus 
FLOMGE (Solid Line) and ISES (Dashed Line) at M = 
0.814, Re = 24.7 x 10 6 , a = 0° 
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Figure 4: RAE 2822 Case 6 Surface Pressures; Com- 
parison of FLOMGE with Johnson-King Model and Ex- 
periment [15] 


There are, however, data cases which stress the credi- 
bility of all currently available 2-D airfoil codes. An ap- 
parently innocuous example is provided by the NACA 
0012 airfoil at zero incidence. This condition removes 
the angle of attack as a “Fudge Factor*. In figures (5) 
and (6) we compare experimental data [15] with results 
computed by FLOMGE and ISES (viscous/inviscid cou- 
pled code developed by Giles and Drela [13]) at two 
different Mach numbers. The solutions computed by 
FLOMGE and ISES agree well with each other at the 
lower Mach number, but the computed shock locations 
are too far back on the airfoil. At the higher Mach num- 
ber, the ISES result is a little better than the FLOMGE 
solution, but again the shocks are too far back. The 
transition point for these calculations was placed at 
3% of chord, but changing the transition point location 
drastically, say to 40% of chord, changes the shock lo- 
cation very little. These poor test/theory comparisons 
are present not only for ISES and FLOMGE but for all 
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Figure 6: NACA 0012, experiment (Symbols) versus 
FLOMGE (Solid Line) and ISES (Dashed Line) at Af = 
0.835, Re = 24.7 x 10 6 , a = 0° 
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Figure 9: Skin Friction Computed on (256 x 64) through 
Figure 7: Pressure Distributions Computed on a Se- (1024 x 256) Grids for RAE 2822 Case 7. 
quence of Grids (64 x 16) through (1024 x 256) for RAE 
2822 Case 7. 
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Figure 8: Pressure Distributions Computed on (256 x 
64) and (1024 x 256) Grids for RAE 2822 Case 7. 


Figure 10: Skin Friction Computed on (256 x 64) 
through (1024 x 256) Grids for RAE 2822 Case 7 (Focus 
Near Leading Edge). 








Navier-Stokes codes we have tried and for all turbulence 
models. 

We realize that one must properly account for wind 
tunnel effects, especially for 2-D flows. However in 
the NACA 0012 test case it would require a Mach 
number reduction of more than 0.02 to produce a rea- 
sonable test/theory comparison. The large number of 
comparisons we have made between experimental data, 
ISES, and available 2-D Navier-Stokes codes suggest 
that there really is something wrong with the Navier- 
Stokes codes and/or the wind tunnel data which must 
be corrected. 

2— D Grid Refinement Studies 

As part of our program to find out if something is 
ailing the 2-D Navier-Stokes codes, we have taken ad- 
vantage of the large memory afforded by our CRAY Y- 
MP to make exhaustive grid convergence studies. We 
have conducted such studies using the central-difference 
scheme of Martinelli and Jameson [2], A grid refine- 
ment study for RAE 2822 case 7 is shown in figures 
(7) through (10). The flow conditions for this calcu- 
lation were taken to be M = 0.73, a = 2.0° and 
Re — 6.5 x 10 6 (based on chord). Transition was set 
at 3% of chord. 

Trailing Edge Glitches 

Glitches in the solution at the trailing edge are quite 
apparent. These glitches are characteristic of JT-NS 
codes for airfoils with a finite trailing-edge angle. The 
glitches do not go away with grid refinement. If any- 
thing, they tend to increase in amplitude. We have not 
found a satisfactory cure for these glitches but they can 
be ameliorated by turning off the fourth order artificial 
dissipation near the trailing edge. 

Martinelli Compromise 

On the 512 by 128 mesh (10), the transition from 
laminar to turbulent flow takes place over about 10 grid 
points. This spreading out of transition is caused by the 
“Martinelli Compromise” in the artificial dissipation, 
which has become common practice in JT-NS codes. 
The “Martinelli Compromise” is introduced to enhance 
convergence on grids with high-aspect ratio cells char- 
acteristic of a Navier-Stokes calculation [2]. Since JT 
codes depend on explicit time marching, the local time 
step they are permitted to use depends on how long it 
takes information to traverse the cell in the short direc- 
tion. For a high aspect ratio cell this does not provide 
time for information to traverse the cell in the long di- 
rection. To ensure convergence, Martinelli dissipates 
the information traveling in the long direction by aug- 
menting the artificial dissipation in this direction. The 


factor by which the artificial dissipation is augmented is 
proportional to the cell aspect ratio raised to the 2/3’s 
power. Thus in the case of a 1000-to-l aspect-ratio cell, 
the artificial dissipation in the long direction is multi- 
plied by 100. 

The “Martinelli Compromise” does serve to improve 
the speed and reliability of convergence, but as can be 
seen clearly at the transition point, the quality of the 
solution is indeed compromised. In some JT codes the 
ill effects of augmenting the artificial dissipation are 
diminished by reducing the 2/3’s power to something 
smaller like 1/2 or even 0.3. The artificial dissipation in 
ARC2D [16] is essentially the same as that present in 
JT codes, except that no compromise is introduced. As 
a result, transition in ARC2D typically takes place over 
3 points. On the other hand it has been our experience 
that ARC2D does not converge as reliably or as fast as 
JT codes. 

Laminar Flow Convergence 

In looking at figures (7) through (10) one notices that, 
as grid density is increased, the airfoil surface pressure 
distribution first begins to lock onto its grid converged 
values (with the exception of the immediate shock re- 
gion), next the turbulent skin friction distribution locks 
in (but not at the shock), and finally (on unacceptably 
fine grids) the laminar skin friction distribution begins 
to lock in. We find it disturbing that a Navier-Stokes 
code would have so much trouble with laminar flow, par- 
ticularly when compared to the resolution requirements 
for accurate solutions in boundary-layer codes. 

For typical airfoils, the boundary layer is laminar for 
only a few percent of chord, and poorly resolved lami- 
nar regions often have little impact on the lift and drag 
calculations. However, we are also concerned with sit- 
uations where laminar flow and transition prediction 
are important; hybrid laminar flow control and high-lift 
devices are two examples. For these flow fields, accu- 
rate prediction of laminar boundary layers on reason- 
able grids is crucial. The behavior at the shock (skin- 
friction reversal only on the 1024 by 256 grid) could also 
have an impact on the pressure drag. 

SAD Laminar Flow Test Case 

The poor performance of methods using scalar artifi- 
cial dissipation (SAD) for high Reynolds number lam- 
inar flows can be demonstrated by considering flow 
over a flat plate at zero incidence. We present re- 
sults for a laminar flat plate at a Reynolds number 
of Re = 500,000 and free stream Mach number of 
Mqo = 0.3. More detailed results for this test case will 
be presented elsewhere. 

Two numerical schemes are employed to solve this 
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Figure 11: Grid Convergence of Velocity Profiles at 

x = 1 for Central Difference with Scalar Artificial Dis- 
sipation 
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Figure 12: Grid Convergence of Velocity Profiles at 

x = 1 for Second-Order Upwind 


Table 1: Central- Difference Boundary Layer Parame- 

ters at x == 1 


grid 

Cj 

( x 1000) 

6* 

(xlOOO) 

e 

(xlOOO) 

C j 

16 x 8 

1.683 

6.240 

3.514 

2.956 

32 x 16 

1.693 

5.343 

2.741 

2.321 

64 x 32 

1.124 

2.977 

1.158 

0.651 

128 x 64 

0.965 

2.564 

1.009 

0.487 

Blasius 

0.9390 

2.434 

0.9390 

0.4409 

Table 2: U] 

pwind Boundary Layer 

Parameters at x — 1 

grid 

Cf 

( x 1000) 

6* 

(xlOOO) 

e 

(xlOOO) 

C/Re$ 

16 x 8 

1.006 

2.435 

1.096 

0.5512 

32 x 16 

0.901 

2.509 

0.994 

0.4481 

64 x 32 

0.919 

2.465 

0.956 

0.4392 

128 x 64 

0.9318 

2.450 

0.9471 

0.4413 

Blasius 

0.9390 

2.434 

0.9390 

0.4409 


flow. The first utilizes Jameson technology — central 
differencing with scalar artificial dissipation — to dis- 
cretize the inviscid fluxes. The second scheme dis- 
cretizes the inviscid fluxes using Roe’s flux-difference 
splitting with second order upwind extrapolation of cell- 
centered states to cell faces [19]. Both schemes discretize 
the viscous fluxes using central differencing. 

Figures (11) and (12) show velocity profiles at x = 1 
unit downstream of the plate leading edge computed 
using the two numerical schemes. The profiles are com- 
puted on a sequence of four grids obtained by deleting 
every other grid line from the finest grid in the typical 
multigrid fashion. The finest grid contains 64 cells nor- 
mal to the plate with the upper boundary at approx- 
imately three boundary layer thicknesses; the grid is 
parabolically stretched away from the plate. The grid is 
also parabolically stretched away from the leading edge 
in x with a grid spacing of approximately Ax = 0.03 at 
x — 1 . 

Figures (11) and (12) show much faster grid conver- 
gence for the profiles computed with the upwind scheme. 
The central-difference results are characterized by an 
overshoot in the velocity near the edge of the bound- 
ary layer and a significant thickening of the boundary 
layer. The two coarsest upwind profiles also show an 
overshoot, but that for the 16-cell grid is no worse than 
the result for the central difference scheme on the 64-cell 



finest grid. 

The disparity in accuracy between the central- 
difference and upwind solutions is further shown in Ta- 
bles 1 and 2, where skin friction, displacement and mo- 
mentum thicknesses are compared with the Blasius pro- 
file parameters at x = 1 . Table 1 shows a quite rapid re- 
duction in errors for the central-difference scheme with 
increased grid density (better than second order), but 
coarse grid errors are enormous compared to the upwind 
scheme results. For the 32-cell grid, the upwind solution 
contains approximately 18 cells within the boundary 
layer and gives 2% errors in the predicted parameters. 
This is consistent with our experience on resolution re- 
quirements for boundary layer solvers. In comparison, 
the central-difference results are still in error by 20% on 
this same grid. 
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Figure 13: x-Momentum Equation Budget at x = 1 for 
Central Difference Scheme (64-cell grid) 

We have identified the culprit for the relatively poor 
performance of the central-difference scheme; it is the 
scalar fourth-difference artificial dissipation in the nor- 
mal direction; it is not related to the “Martinelli Com- 
promise”. Specifically, contamination results from ex- 
cessive dissipation normal to the boundary layer in the 
x-momentum equation. This occurs because the artifi- 
cial dissipation is scaled by the flux Jacobian spectral 
radius \v\ -h c, whereas a properly formulated matrix 
dissipation or upwind scheme (e.g. Roe’s flux-splitting) 
scales the normal dissipation by \v\. It is easily shown 
that with the \v\ + c scaling, the normal artificial dissi- 
pation in the x-momentum equation is proportional to 
(A y/6) 3 VRe/M based on edge conditions. Therefore, 
as the Reynolds number is increased, more grid. reso- 
lution (i.e. more grid points across the boundary layer 
thickness 5) is required to achieve a given level of accu- 



Figure 14: x-Momentum Equation Budget at x = 1 

for Central Difference Scheme (64-cell grid) Blow-up of 
Boundary-Layer Edge 
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Figure 15: x-Momentum Equation Budget at x = 1 for 
Upwind Scheme (64-cell grid) 



racy. This also explains why similar poor performance 
of scalar dissipation methods is not seen in low Reynolds 
number flows. 

To illustrate the contamination, the budget for the 
x-momentum equation for the profile of cells at x = 1 
is plotted in Figures (13), (14) and (15) for the central- 
difference and upwind schemes on the finest grid. In the 
figures, EULI, VISI and DSPI represent the difference 
in inviscid, viscous, and artificial dissipation fluxes, re- 
spectively, through the vertical faces of each cell (i.e. 
streamwise fluxes). EULJ, VISJ and DSPJ represent 
the analogous flux differences through horizontal faces 
(i.e. normal fluxes). For the upwind scheme, DSPI and 
DSPJ are taken as the difference of the split-fluxes and 
the face-averaged fluxes; hence, EULI and EULJ are 
consistently defined between the upwind and central- 
difference schemes. 

Figures (13) and (14) reveal that the normal artifi- 
cial dissipation (DSPJ) is large everywhere in the pro- 
file, even outside the boundary layer. Near the wall, 
the momentum balance is completely nonphysical with 
artificial dissipation (DSPJ) balancing viscous diffusion 
(VISJ). The budget for the upwind scheme is more phys- 
ical; artificial dissipation is small everywhere, and the 
dominant terms are EULI, EULJ, and VISJ. 

Some previous researchers have introduced ad hoc 
scaling reductions of the artificial dissipation through 
the boundary layer as an attempt to eliminate contami- 
nation. We have also applied some of these “fixes” with 
disappointing results, and know of no ad hoc scaling 
that will reduce the artificial dissipation across the en- 
tire profile to a point where the results are comparable 
to a properly formulated upwind scheme. 

We wish to emphasize that these problems with 
high Reynolds number laminar flows are not inherent 
to central-difference schemes, but to central-difference 
schemes that use scalar artificial dissipation (CDS AD). 
This leads us to conclude that any scheme using scalar 
artificial dissipation or any scheme that is highly dissi- 
pative for low Mach number flows (e.g. van Leer’s flux- 
splitting, see Ref. [5]), should be suspect for calculating 
laminar flows. 

Our current research is directed towards improving 
the convergence rates of upwind schemes to steady- 
state. It is well known that reducing the spatial dissi- 
pation in a scheme usually results in slower convergence 
to steady-state. 

2-D High-Lift Configurations 

High-lift flow provides a significant challenge to CFD 
technology. For instance, the CFD code must have the 
ability to accurately predict the laminar boundary-layer 
profile ahead of the transition point so that a transi- 


tion prediction method can be applied. The confluent 
boundary layer on the main element and the separated 
flows in the cove and on the flap must be modeled. 
There are free-shear layers in many parts of the flow- 
field where the spatial length scales of the flow char- 
acteristics are non-isotropic. The free-shear flows in- 
teract with the boundary layer on the flap to some- 
times cause dramatic and unexpected flow behavior 
(e.g. Reynolds number reversal effects described in [8]). 
Simple boundary-layer approximations may not be ad- 
equate for such complex flows. Navier-Stokes methods 
seem to be the natural choice, but even here turbulence 
models remain a major issue. 

We have written a code called A610 described in [17] 
that uses viscous/ inviscid coupling to calculate flows 
around multi-element airfoil configurations. We will 
compare A610 with the Mavriplis unstructured grid NS 
code [18] for a Douglas 3 element configuration tested 
at LTPT. Comparisons between experiment, A610 and 
the Mavriplis code for 8, 20, and 23 degrees angle of 
attack are shown figures (16), (17), and (18). In order 
to run with A610, the coves on the lower surfaces of the 
leading edge slat and near the rear of the main element 
had to be smoothed. The effects of this smoothing are 
particularly noticeable in the A610 results at 8°. At all 
angles of attack A610 seems to predict C p peaks that 
are a little too high. The overall test/theory compar- 
isons seem to favor A610 at 20° and the Mavriplis code 
at 23°. At 8° A610 properly predicts separation for the 
trailing edge of the flap while the Mavriplis code does 
not. The inability to predict this flow separation seems 
to be a failing of the Chimera based Navier-Stokes codes 
as well. 


Practical CFD assessment for High-Lift 

Given these results there does not appear to be any 
strong reason for us to favor the Navier-Stokes code. 
All the more so since we know that being a NS-CDSAD 
code, the Mavriplis code is not able to properly cal- 
culate the laminar portions of the boundary layer and 
thus can not give us a transition prediction capability. 
The importance of transition, shown in figure (19), is 
computed using A610. When the Navier-Stokes codes 
come closer to achieving their theoretical potential, we 
will use them in earnest. 

Also, while the preliminary capability in 2-D is being 
developed by many researchers, we badly need a 3-D 
code. In three dimensions, high-lift flow can be even 
more complex than in two dimensions. The edge vor- 
tices, gap flows, and embedded longitudinal vortices in 
the boundary layer all have strong effects on the overall 
performance of the high-lift system. 
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Figure 17: “High lift Olympics” 20° angle of attack. 
Experiment (Symbols) vs. Mavriplis (Lines) and A610 
(Dashed-Lines) 


Figure 19: High lift, importance of transition 


Nacelle-Flow Analysis 


Nacelle analysis and design is an integral part of the 
airplane design process. In advanced aircraft, propul- 
sion systems are closely coupled with the airframe, and 
proper engine installation is essential in order to im- 
prove the overall performance of the aircraft. Inviscid 
methods (panel and full-potential) have been very use- 
ful, but viscous effects are also of interest especially for 
off-design conditions on large twin-engine airplanes. 

It is relatively easy to analyze an isolated flow- 
through nacelle using a code written to treat wings. We 
have adapted the TLNS3D Navier-Stokes code. The na- 
celle is treated as a ring wing, with periodic boundary 
conditions. To simulate a powered nacelle one can ei- 
ther specify inlet and exhaust boundary conditions, or 
use a center body with variable geometry to control the 
mass flux through in the engine. At cruise condition, 
the Navier-Stokes code provides accurate results, simi- 
lar to that of wing or wing/body analysis. Problems are 
encountered in nacelle analysis with low-speed takeoff 
conditions, and with high-speed engine-out conditions. 

At takeoff, the effective angle of attack for the na- 
celle is high. The flow is highly three-dimensional, and 
a laminar separation bubble may form at the nacelle 
highlight region. The marginal accuracy of the avail- 
able Navier-Stokes codes in the laminar flow region was 
mentioned above. In the long run we must arrive at a 
reliable 3-D boundary-layer transition-prediction capa- 
bility, as well as a plausible behavior in the transition 
region, before we can capture the laminar separation 
bubble. This bubble has dramatic effects on the overall 
flow field. Figure (20) compares the results of nacelle 
analysis, first treating the flow as fully turbulent (tur- 
bulence model active in the whole domain), and then 
assuming transition at 5% from the leading edge (turbu- 
lence model active only downstream of that line). The 
results are drastically different, and neither agree well 
with experiment. The flow pattern with transition at 
5% is however similar to the experimental pattern. 

At a high-speed, engine-out condition the large 
amount of spillage around the nacelle results in a strong 
shock on the exterior surface of the fan cowl, which may 
cause severe shock-induced separation. Present Navier- 
Stokes technology is capable of handling mild shock- 
induced separation. However, none of the turbulence 
models tested gives reliable solution for strong shock- 
induced separation. 

In summary, attempts at nacelle analysis and engine- 
airframe integration by Navier-Stokes solutions raise the 
same issues as wing design. These are: gridding diffi- 
culties when other components are included; numeri- 
cal accuracy particularly in the boundary layers; and 
turbulence-modeling accuracy particularly at shock in- 
teractions. In addition, because of lower Reynolds num- 


Erfects of trip location on nacelle lip separation 
(High alpha, low Reynolds No.) 



Figure 20: Nacelle with transition at 5% versus all tur- 
bulent flow 


bers and extreme velocity peaks at the lips, laminar re- 
gions may exist in the boundary layers and exert much 
control over the global flow field. In the long run we 
need a reliable and, as much as possible, automatic 3- 
D boundary-layer transition-prediction capability. For 
this, two key ingredients are — presumably — a stabil- 
ity analysis with sufficient robustness and generality 
to handle steep three-dimensional pressure gradients, 
and accurate velocity profiles directly out of the Navier- 
Stokes solver. Neither ingredient is at hand. The turbu- 
lence models also need improvement to handle moder- 
ate or massive separation, whether encountered at low- 
speed takeoff conditions or at high-speed, engine-out 
conditions. 


Conclusion 

The 3-D Wing/Body calculations show that Navier- 
Stokes codes hold much promise. However, our 
test/ theory comparisons in 2-D and for nacelles, as well 
as our detailed 2-D grid refinement studies, are sober- 
ing. It is apparent that much work remains to be done 
in numerics and physical modeling of transition and tur- 
bulence before we can say that we have an “Industrial- 
Strength” Navier-Stokes code in hand. 
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